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Abstract 

We consider an XXZ spin- 1/2 chain in the presence of several types of disorder that do not break 
the XY symmetry of the system. We calculate the complete asymptotic form of the spin-correlation 
functions at zero temperature at the transition between liquid and disordered phase that occurs for 
a special value of anisotropy in the limit of small disorder. Apart from a universal power law decay 
of correlations, we find additional logarithmic corrections due to marginally irrelevant operator of 
disorder. 

PACS numbers: 75.10.Pq,71.10.Pm 
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I. INTRODUCTION 



A spin-1/2 XXZ spin chain is an exactly solvable model [l| and has the Hamiltonian 

#xxz = JJ2 [i S t S i+i + S i S i + i) + *S' SUi] > (!) 

i 

where the spin can be taken as S = er/2, where a are the Pauli matrices. They satisfy 
the commutation relation [S" 1 , S%] = ie a b c S^5jk- XXZ model exhibits rich physics when 
varying the anisotropy parameter A. For A < — 1 the ground state has Ising-type long- 
range ferromagnetic order and finite excitation gap in the low-energy spectrum. For A > 1 
the ground state is again gapped and posses antiferromagnetic order. For — 1 < A < 1 the 
model is gapless and has a linear spectrum of low-energy excitations. For — 1 < A < 1, the 
staggered part[2[ of its asymptotic correlation functions read [l|, 0] 

W> ~ (-i) j br 1/2i s (2) 

(s»s]) ~ (-iy\j\- 2K , (3) 

where the parameter K is given as 

A = - cos ^- (1/2 < K < oo). (4) 
2K 

For A = 1 the model ([1]) becomes isotropic and thus the correlation functions ([2]) and 
become equal. However they acquire an additional logarithmic factor 



• 1oe 1/2 I ? I 

(sssf) = (s z s]) ~ (5) 

Its origin is well understood. The low-energy properties of (pQ) correspond to a (1 + 1)- 
dimensional sine-Gordon (SG) field theory. For A = 1 the corresponding SG model is 
exactly at criticality where its cosine term becomes marginally irrelevant and produces the 
above logarithmic corrections. Numerical works 0] have confirmed its presence in (JSJ). 

An insight into log- correct ions of the Heisenberg chain (XXZ for A = 1) comes from 
the fact that the SG model is equivalent to the classical two-dimensional XY model that 
exhibits a Berezinskii-Kosterlitz-Thouless (BKT) transition [l^, EH- In his early work 12 



Kosterlitz showed that the spin-spin correlation function at the transition contains loga- 
rithmic corrections of the form log 1 / 8 |i|/|i| ly/4 (there the spins are classical two-dimensional 
vectors). The ratio of powers of log \j\ and \j\ is 1/2, as in (J5J). That is easily understood 
since the power law exponent fixes the exponent of the logarithmic part for the same type 



of correlation functions, as in the case for the two examples. Having in mind the work [12 
it is surprising that the logarithmic corrections for correlations of the Heisenberg chain were 
rediscovered relatively late. 

In the present paper we consider several types of disordered XXZ chains. The form of 
disorder is such that is preserves the XY symmetry of the pure model (PQ), see below. It turns 
out that all the disordered models we consider here have a BKT transition at A = —1/2 



in the limit of vanishing disorder strength [13l . Il4j , and the region — 1 < A < — 1/2 (where 



the disorder is irrelevant) is a remnant of the region — 1 < A < 1 of the pure model, 
and therefore has quasi-long-range order and correlations (T5]) and ([3]), for weak disorder 
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[l3| | . We will calculate the spin-correlation functions at the transition. As in the pure case, 
they acquire again multiplicative logarithmic corrections, however the universal powers of 
logarithms are now different for (J2]) and fl3]) . That comes from the fact that the corresponding 
field theory is given by a disordered SG model. 

The article is organized as follows. In Section II we introduce the different types of 
disorder and describe the low-energy properties of the system in terms of a 1+1-dimensional 
field theory. In Section III we compute the correlation functions of interest. Section IV 
contains conclusions. In Appendix we give the derivation of the scaling equations for the 
model. 

II. DISORDERED HAMILTON! AN AND LOW-ENERGY DESCRIPTION 

A. Different types of disorder 

We are interested in different types of disorder that preserve XY symmetry of ([I]): (i) a 
random transverse magnetic field along z direction, 

H ZF = KSt, (6) 

i 

(ii) random z— exchange, 

i 

and (iii) random planar exchange interaction, 

i 

Note that h*, 5 J* and 5J* y are random variables uncorrelated from site to site, each of them 
having the zero mean and the second moment determined by D: 

hfh* = 5J*5J* = SJ- y 5J- y = DSij. (9) 

The overbar denotes an average over different disorder configurations. In the following we 
consider the Hamiltonian 

H = H xxz + H D , (10) 
where H D denote any one of the terms (jSD, dZj), or (JSJ). 



B. XXZ chain 



First we start by considering the XXZ chain, which treatment is now standard [la, [16 
One translates the discrete model (CQ) into a continuous field theory that captures long- 
distance physics. Using the Jordan- Wigner transformation one first translates A = part 
of (OP) into noninteracting fermions. The remaining z— part essentially represents the in- 
teraction between fermions, since the spin operator S z becomes the density of fermions. 
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After employing the bosonization representation one ends up with the continuum form of 



the action [la, 16 



S = ^JdxdT 

S g = —g J dxdr cos (4<£>), 
where g ~ J A, and the spin-wave velocity is 



vK{d x 6f + -{d xV f + 2i(d x 6)(d T <p) 



v = Jci- 



Iy 



2K- 1 



sin 7T 



2.RT- 1 
2K 



(11) 
(12) 

(13) 
(14) 



while the Luttinger parameter K is connected to the anisotropy via (j3J). The lattice constant 
is denoted by a. The spin operators are connected to the bosonic fields (p(x,r),9(x,r) 
through the following relations [Ts], [l6[ : 



S 



imkp Xj +im,(p(xj )—i8(xj ) 



, . ., V B 

m even 

S z I , 

-V(p(xj) + } C m e 

IT < * 



3 

a 



7T 



^imkpXj+irnip(xj ) 

m^0,even 

where Xj = ja and fc^ = it /{2a). Also, _B m and C m are some nonuniversal constants. 



(15) 
(16) 



C. Disordered chain and scaling equations 

Next we include disorder. We start with the random field case fl6]) first. The corresponding 
action is 

S Z¥ = J dxdr ^-~rj(x)d x (p + 2C 2 t(x) cos(2^)| + . . . (17) 

where 



7](x)r](x') = D5(x-x'), (18) 



£(x)£(x') = D5(x-x'). (19) 

Here the fields T)(x) and £(x) are proportional to the uniform and alternating parts of the 
random field h z . Higher order terms denoted by ... in (fl~7l) correspond to less relevant higher 
harmonics from ([TBI) . After performing the disorder average using the replica method, we 
find the action that corresponds to H xxz + H ZF : 

s = J2 s o(<Pa,e*) + s f + s b . (20) 

a 

The first term is given by f|T2|) . The off-diagonal quadratic part that comes due to disorder 
is 

5/ = - 2~2' D/ S / ^TdT / [d x (p a {x,T)][d x (pp(x,T / )], (21) 
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while the anharmonic part of the action reads as 



S b = -C\D h / dxdrdr cos[2<^ a (x, r) - 2cp (3 (x, r')]. (22) 

a/3 J 

Greek letters a, /3 — 1, . . . ,n denote replica indices, where the number of replicas n is sent 
to zero at the end of calculations. Note that only the most relevant terms are written, 
since they determine the critical behavior. The term f fl3|) turns out to be unimportant, 
since it is irrelevant close to the transition 13||, occurring for K = 3/2 (A = —1/2). The 



term analogous to ( 12"2"|) but having plus sign in the argument of cosine is also suppressed 
in our calculation. The bare parameters and Df are the same and equal to D in the 
unrenormalized model. However, they renormalize differently, hence we distinguish them. 
The model ( 120]) is well known and describes disordered bosons and fermions in one di- 



mension 



,,,,, ,, L [14jL an d via Jordan- Wigner transformation straightforwardly extends to disordered 
spin chains [13j. It has the following scaling renormalization group equations 14|: 



di 
d5 R 



-2@ R 5 R + (D(@ R ), (23) 
-9@ R + 0(@ R 5 R ), (24) 



di 
-( 

dV~ Jn/ di \K R 

where i denotes the scale and the subscript index R denotes renormalized coupling constants 
that flow and that uniquely correspond to the bare ones. The dimensionless disorder strength 
is introduced as 3) = 7rCf a 3 _Db/t> 2 , while 5 = K — 3/2 denotes the distance from the 
transition, occurring at K = 3/2 or equivalently A = —1/2, see (jl]). For A < — 1/2 weak 
disorder is irrelevant due to quantum fluctuations and at large scales one expects to have 
the properties of the XXZ chain. The phase for A > — 1/2 is disorder-dominated. Due to 
perturbative nature of the scaling equations, the above picture is valid for sufficiently weak 
disorder. Strong disorder even for A < —1/2 is expected to dominate the physics. 
The solution of the flow equations (123]) and ff24l) is 



5 2 R - 9@ R = C, (26) 

where C is an arbitrary constant, and C < (C > 0) marks the disorder dominated (quasi- 
long-range ordered XXZ) phase. For C = the system is at the transition, where the 
solution of (JMD is 

8*W = ^ (27) 



III. CORRELATION FUNCTIONS 



In this section we calculate the spin-correlation functions {SfS~ + j) and (S?S? + j) at the 
transition between liquid and disordered phase, occurring for C = in (T26]) . We focus below 
on their staggered part and moreover only on the most relevant contribution. We briefly 
discuss the uniform part of correlations at the end. 
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Using (fT5l) . and (|T6|) we have: 



(StS~ +j } ~ (-l)^^! ^' )- ^. )]) = (-l) J 'ii! + _O,0), 



(SfS^) ~ (-iya 2 C 2 2 <e i2 ^+^°)-^' )] +h.c.) 
= (-l) J a^[4(ja,0) + k], 

The average is with respect to the full action Sq + S ZF , Eqs. f|T2|) and ([[ 



(28) 
(29) 



A. R zz (x,t = 0) correlation function 



The Callan-Symanzik equation determines the critical correlation function as [17 

R ab (r, 0; 5, 0) = e^ r d ^W# afe (a, 0; @ R {t r )), 
where l T = In (r/a) ^> 1. Here 

9 In Z ah 



lab 



d In a 



(30) 



(31) 



where Z a ;, is the prefactor that multiplicatively removes the a-dependent divergences in 
R a b- At large scales R a b(o>, 0; S R (£ r ), S> R {£ r )) becomes a constant since the parameters 
8 R (£ r ), @ R {£ r ) reach their fixed point vales, that is zero in our case. Note that everywhere 
below in the text we use simplified notation R a b(r, 0) = R a b(r, 0;5,@). 

We first calculate R zz (r,r = 0) (r > 0) using the perturbation theory in the disorder 
strength [181 ] . Using the replica action we do the disorder average and should evaluate 



R zz (r,0) = ii m (e*%7M)-^(o,o)]^ 



n-K) 



where S is given by (TSUI) . After evaluating Gaussian integrals one obtains 



where W\ is 



Wi(r) 



#z*(r,0) = e" 



C 2 a 3 D b 



3+25 



[1 + Wi(r) + 0(^<J)] 



(32) 



(33) 



dxdydy' 



[a 2 + {y-y' 



,2l3/2 



Ai + a 2 A 4 + a 2 
A 2 + a 2 A 3 + a 2 



3/2 



9 



^ln 2 



2+C 2 ^hl— . 

4 er cr 



(34) 



and C2 being some number. We have introduced the abbreviations 

A x = x 2 + y 2 , A 3 = (x 



r) 2 + y 2 , 



A 2 = x 2 + y' 2 , A A = {x- r f + y' 2 . 



(35) 



W\ contains logarithmic divergencies for a — > 0. In order to calculate the long-distance 
behavior of R zz we should replace the bare parameters by the renormalized ones, which are 
connected via the relations: 



5 = 8 R + 9( Cl - \)® R /2 + 0(& R 6 R ), 
@ = @ R + 0(@ R 5 R ). 



(36) 
(37) 
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Here A = log (c 2 m 2 a 2 ) and c is the constant. The details of calculations are presented in 
Appendix, see Eqs. (IA7j) and (1A8I) . Using (1361) and (|37|) one ends up with 



i2zz(r,0) =e" 



2K 2 -D^ / r \-3-2<5 H 



1 - \®r\ 2 + ^ (9 Cl - 2c 2 ) ^ fl A 



x {1 + [(2c2 - 9ci) ln(cror) + 91n 2 (cmr)] + 0(#ji<Jr)} 
so R zz is renormalizable by the multiplicative factor 



[ma) 



-3-25 



R 



\ 9 

1 + 4^ 



A 



(9ci - 2c 2 ) ^ R A 



(38) 



(39) 



Importantly all r-dependent divergencies for a — > have canceled once one uses the renor- 
malized parameters. Then one finds the anomalous dimension function 



d In a 



-3-2S R + 0(9 R ,8 R ). 



Finally, using Eqs. (127]). (130]) . and ( HP]) we find the large scale behavior 



Rzz(r,0) 



~ e 



*)V 2 ( T - 
r J \a 



(40) 



(41) 



where we used the fact that Dj and K/v are not renormalized, see f )25p . We have found a 
logarithmic correction to the correlation function at criticality. Note that higher order terms 
in (j4"0"]) do not give important contribution since after integration in fl30|) they saturate into 
a constant for large x. This can be seen as a change of the number under the logarithm that 
we neglected. The same type of contribution comes from (1271) at short length scales, while 
the universal large scale part l/£ determines the logs we find in (j4ip and below in (I4"9l . 



B. (r, r = 0) correlation function 



Next we calculate the other correlation function 



J R + _(r,0) = lim(e^( r '°)^( ^) 5 . 



n->0 



(42) 



Since it contains the conjugate field to (p, the disorder term proportional to Df does not 
appear in any order of perturbation theory in One way to see that is to perform a 

gauge transformation by absorbing Df term in the new tp fields. Under such transformation 
(|22|) is invariant. After performing the perturbation theory one obtains 



Cl \ 3+25 
V 



[I + W({t) + 0{@S)] 



(43) 



where 



W?(r) 



Cia 3 D b 



1 

4'' 



dxdydy' 



[a 2 + (y-y' 



,21 3/2 



^-.2 ^,2 

'In 2 — + c 3 ^ln— . 

a z a z 



MM i (y 



y'fr 2 



A 2 A 3 2 V / A 1 A 2 A 3 A 4 



- 1 



(44) 
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Here c 3 is some constant. Note that A±, . . . A4 have been defined in (1351) . Then one gets 



R+-{r,0) 



3+26 



Plugging in the renormalized parameters 

1 



'ln 2 - + c 3 ^ln— + 0($>6) 
and ( 157)) one gets 



(45) 



#2 r 



3 + 2Sr 



1 + ^A 2 



-(2c 3 + Cl )^A 



[l - ^ln 2 (cmr) + (2c 3 + d) 3> R \n(cmr) + C(^r(5r 



(46) 



One should notice that all r-dependent divergencies in a — > limit of the perturbation 
theory in (j4"3"]) cancel once one expresses the bare coupling constants by the renormalized 
ones, see Eq. (1461) . We have the combination of the ln 2 (r/a) term from the right hand side 
of ()45p that combines with — ^Aln(r/a) that arises after expanding the power in (143)) with 
the exponent 1/(3 + 28) after using the renormalized quantities. 

The previous expression (T4"6")) contains divergence when a — > that could be removed 
multiplicatively by 



J R4 



[ma) 



3+2S R 



(47) 



so that the renormalized correlation function Z R+ _ 
anomalous dimension function then reads 

d\nZ R+ _ 1 2 



R + - (r, 0) has no dependence of a. The 



1+- 



5 R + 0(® R ,6 2 R ) 



(48) 



din a 3 9 

where we have used the scaling equations ( 123)) . and ()24j) . Then, using Eqs. (1271) . (15D1) . and 
151) we find the large scale behavior 



«-M>~(;) l/V/s G 



(49) 



From ( 1281) . ( 1291) . (14T1) . and (1491 we obtain the equal time spin-spin correlation functions 
for H xxz + if zp at the transition (near A = —1/2 and for weak disorder): 



X 

i+j 



) 



2(SfS; 

(~1) J 
|j| 3 ln 2 \j\ 



-l) J ln 2/9 |j| 
|j| 1/3 



(50) 
(51) 



using the exact value at the transition v/K = \/3Ja/4. We notice that the uniform non- 
staggered part [15j for (S^~S~ + j) is subleading at large distances (~ |j| -10 ^ 3 ) therefore neg- 
ligible with respect to the staggered part (15T)1) . However, the uniform part for (SfSf + j) is 
~ |j|~ 2 and therefore dominates at large distances. Finally we notice that inside the logs of 
the last two formulas in general there are some nonuniversal disorder dependent numbers 
that come from small scales of the exponent in the Callan-Symanzik equation (130)) . that we 
neglected at large distances but that may be important for numerical analysis. The simi- 
larity of effective models for disordered spin chains and disordered bosons (and fermions) 
in one dimension directly determines the form for the single particle and density-density 
correlation functions for the latter, that respectively have a form similar to (150)) and ( 15~T|) 



19 
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C. Random exchange models 



Next we discuss disordered Hamiltonians given by (J7J) and (JSJ). At first glance they appear 
quite different from the random field case ([6]); however due to quantum fluctuations the 
critical phase (for —1 < A < —1/2 and small disorder) has the same properties for all three 
cases, (j7j), flS]), and (jSJ) [13]. That occurs because the disorder is irrelevant there. On the 
contrary, properties of disordered phases are different, since they are sensitive to the form of 
disorder. It turns out that critical properties at the transition from a phase with quasi-long- 
range order to a disorder-dominated phase are the same for all three cases [l3[. However, 
the corresponding replicated action for © and (jSJ) is given by S = J2 a So(<Pa, a ) + Sb- We 
see that the term (j2ip . or equivalently the first term in ffTTl) of the unreplicated action, is 
missing for the two random exchange models 13, 2(j. That has very important consequences 



for the (S?Sf + j) correlation function, since it does not get suppressed by the disorder as in 
f[5"T]) . Therefore, one can repeat the above analysis and obtain the correlation functions at 
the transition. The first one will have the form as in (|50|) . while the other one is now changed 
and reads 

Ijfln 2 



(S!St +j ) ~ (52) 



The absence of disorder in (J52l) makes it universal at the transition, which one may detect 
in numerical simulations. 



IV. CONCLUSIONS 

In this paper we have calculated exact asymptotic form of the spin-correlation functions 
for several types of disordered XXZ spin chains at the transition between the critical (dis- 
order irrelevant) and disordered phases at zero temperature. We found the corresponding 
logarithmic corrections, given by (|50|) and (|51"|) for the random field case ([6]), and by (|50|) and 
(152]) for the random exchange cases <^$) and (jHJ). They appear due to marginally irrelevant 
operator of disorder, quite similar to the Heisenberg spin-chain case. While the prefactors 
in the correlation functions for the Heisenberg case are exactly known [2l|, 22| . its precise 



determination in ( [50]) -( j52]) for the disordered case is a formidable task. However, if the 
disorder is weak one may expect the prefactors to have similar values as the ones of the pure 
XXZ case at A = —1/2 that are exactly known [23[. Finally we mention that higher order 
scaling equations [l9| of the model fl20|) would produce additive subleading terms in (150]) - 
(j52p . and therefore they are inessential. The presence of log-corrections at the transition 
may be important for numerical studies for precise determination of the phase boundaries, 
for example. They will also play a role in finite temperature quantities such as the spin 
susceptibilities and the NMR relaxation rate. 
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Appendix A: Effective action and scaling equations 



In this appendix we derive the scaling renormalization group equations for the model 
fl20|) . To treat the backward scattering part of the disorder, we use a field theoretic approach 



24j . |25| to obtain the effective action of the model T((p). It reads 24 



-r(¥>) 



Vxe 



sr 

Sip(x) 



(Al) 



This general equation can be solved perturbatively, order by order with respect to the small 
parameter D b [l9(. Up to an additive constant, the effective action reads 



S + S f + -(S b ( V + X )) x + O(Dl 



(A2) 



where (• • • ) x denotes an average with respect to the field \ using the quadratic part So + <Sf- 
At the lowest order we easily obtain 



Ti = - 



/ dxdrdr' 

a/3 J 



3 4G(0,r-r') _ 1 



>a/3 



+ l|cos[2^ Q (x,r) -2cpp(x,r% (A3) 



where SB = ClD b e~ AG ^\ while 



K 



G(x,t) = —Kq(i7iVx 2 + v 2 r 2 + a 2 ). 



(A4) 



One should notice that only a part of the full correlation function that is diagonal in replica 



indices, (1A4|) . enters the result. This is due to the fact that the off-diagonal part of the full 
propagator does not depend on imaginary time. We introduced m as an infrared cutoff by 
adding a term oc m 2 ((p a ) 2 into Sq, in order to calculate the effective action. 

The obtained perturbative expansion of the effective action contains all the information 
about critical properties of our model at the first order in D b . In order to derive the scaling 
equations one should find the divergent terms in T in the limit a — > 0. After expanding the 
operators one obtains the effective action 



T=J2 [ dxdr{^[(9,^) 2 + m 2 (^) 2 ] + 
- tt4 ^ / dxdrdT'[d x ip a (x,T)][d x (pp(x,r')] 



2tcKv 



(d T ip c 



(A5) 



The coefficient a\ arises from the expansion of operators in (IA30 and reads d\ = 
J drr 2 [e 4G (°' T ^ — l] . Introducing the small parameter 5 = K — 3/2 which measures the 
distance from the critical point and using flA4j) we get 



a i 



-A + ci + Q(5) 
(cmv) 3 



(A6) 
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where A = In c 2 m 2 a 2 , c = e 7B /2, 7^ is the Euler constant, and c\ is a constant. 

The effective action ( 1A5j) has divergencies when ma — > contained in A. In order to 
remove them we introduce a set of renormalized coupling constants (denoted by the subscript 
R ) by & = v = Zv R , 5 = Z(3/2 + S R ) - 3/2, K = Z(3/2 + S R ), m R = m where 

Z b =1 - 5 R X, (A7) 
Z =1 - 3@ R \ + 3a@ R . (A8) 

Differentiating the bare coupling constants with respect to the scale t = — In m R one obtains 
the scaling equations fr23|) -( |25|) . 
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